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^ Abstract 
O 

We study the spectrum of scalar I = quasi-normal frequencies of Anti-de-Sitter 
hairy black holes in four and five dimensional designer gravity theories of the Einstein- 
scalar type, arising as consistent truncations of A/" = 8 gauged supergravity. In the 
— I dual field theory, such hairy black holes represent thermal states in which the operator 

corresponding to the bulk scalar field is condensed, due to the multi-trace deformation 
I associated with non-standard boundary conditions. We show that, in a particular class 

of models, the effective potential describing the vacua of the deformed dual theory can 
be identified, at large values of the condensate, with the deformation plus the conformal 
coupling of the condensate to the curvature of the boundary geometry. In this limit, we 
show that the least damped quasi-normal frequency of the corresponding hairy black 
^ holes can be accurately predicted by the curvature of the effective potential describing 

^y-^ the field theory at finite entropy. 

OO 

o 



^ 1 Introduction 

> 

• ^ The last stage in the approach to a black hole geometry consists of a superposition of 

^ damped oscillations, whose complex frequencies are called quasi-normal frequencies. Fol- 

lowing the seminal work [1], increasing efforts have been devoted to the computation of 
quasi-normal frequencies of "black" geometries with AdS boundary conditions, motivated 
by the AdS/CFT correspondence. Indeed, when a dual field theory description is available, 
we know that these frequencies govern the approach to thermal equilibrium, appearing as 
poles of the finite temperature retarded correlators of the dual operator [2, 3]. 

In this paper, we study the quasi-normal frequencies of four and five dimensional black 
holes dressed with a single neutral scalar field, motivated by consistent truncations of = 8 
gauged supergravity. The scalar field mass squared belongs to a particular range in which 
both asymptotic modes at the AdS boundary are normalizable, and one is free to impose 
general boundary conditions that mix their coefficients, denoted (a, f3) in the following. 
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The resulting Einstein-scalar theories are called designer gravity theories, because their 
properties depend on the free choice of a function that determines the boundary conditions. 
In the dual field theory, such function appears as a deformation of the action: this freedom 
made it possible to engineer holographic setups opening new perspectives in traditionally 
hard research topics, like the formation and possible resolution of cosmological singularities 
[4, 5], and strongly coupled field theories at quantum critical points [6]. 

Effective potentials built from supergravity data are an important tool in the analysis 
of these setups [7]. They are computed from the coefficients (a,/3) of the asymptotic modes 
of the scalar field that characterize families of hairy black hole solutions, and include the 
deformation which corresponds to the boundary conditions. Starting from the family of 
hairy black holes of a given radius Vh, one can build a function Vr-f^ict)^, whose extrema 
a* correspond to the hairy black hole solutions of that radius that fulfill the boundary 
conditions. Moreover, Vr^ia*) gives the total mass of such solutions. In the dual field 
theory, a represents the expectation value of the operator that condensates in the hairy 
black hole states, while the mass and the horizon area are naturally identified with the total 
energy and the entropy at equilibrium. Therefore, Vr^ is interpreted as the effective field 
theory potential describing the thermodynamics at fixed entropy: the equilibrium states are 
extrema of the total energy. An analogous construction starting from the family of hairy 
black holes of given temperature T leads to a function Vt that is naturally interpreted as 
the effective field theory potential in the canonical ensemble at temperature T. Indeed, the 
extrema a* correspond to the black hole solutions of temperature T that fulfill the chosen 
boundary conditions, and VT(a*) gives the free energy of such states. 

At the level of computing the energy/free energy of the equilibrium states, these effective 
potentials encode all the information about the strong interactions that characterize the 
dual field theory: they appear to provide a genuinely non-perturbative picture of arising 
from supergravity data. Previous studies qualitatively showed that they also give a correct 
intuition for the stability/instability of the equilibrium states under perturbations of the 
condensate, depending on the boundary conditions [8]. Our work aimed at generalizing this 
observation, and verifying whether precise information about the quasi-normal frequencies 
of the hairy solutions can be extracted from the form of the effective potentials themselves. 

In Section 2, we introduce the general bulk setup, consisting of a single neutral scalar 
field with negative mass squared coupled to Einstein gravity, and review the existence and 
properties of AdS static hairy solutions, as well as of the effective potentials. We also recall 
the definition of quasi-normal frequencies, briefly describing the numerical procedure we 
have adopted. In Section 3, we present numerical results for five dimensional hairy black 
holes, choosing regularized boundary conditions of the type considered in [9]. In Section 4, 
we generalize our results to a four dimensional model that exhibits similar properties and, 
based on a further four dimensional model, we discuss the limitations of our approach to 
quasi-normal frequencies based on the effective potentials. 

^In the degenerate case = m%p, V,,, is more conveniently expressed as a function of the other 
coefficient, denoted by /3 in the following 
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2 Einstein— scalar setup 



Spherically symmetric, static solutions 

We start from the general Einstein-scalar action: 

1 ^ -id^ I— 1^^5..^2 



The potential V is taken to have a negative local maximum in (/3 = 0, corresponding to a 
unit radius AdS solution: 

Denoting vn?^p the Breitnerlohner-Freedman bound, rrP'^p = —{d — 1)^/4, we will consider 
models with rri^p < m? < m^p + 1, for which both the asymptotic modes of the scalar field 

^BF 

coordinates the AdS metric reads: 



are normalizable. The degenerate case = rn'^p must be considered separately. In global 



ds^ = -(1 + r^)dt^ + ^ + r'dnd^2 

1 + 

and the scalar field decays near the AdS boundary r = oo as: 

ra^p < m? < m^p + 1 : ip = a r~'^^ + (3 r~'^- 
= rri^p : if = a r^^ log r + /? r 

d-l 



(2) 



In this range of masses, basic physical requirements such as the finiteness of the total energy 
do not require the slower decaying mode to vanish, a = 0. Therefore, we are free to consider 
more general boundary conditions, that are conveniently parametrized as: 



m%p Km^ <m\p + l: f3 = W,a 

2 2 

m = m^p a = —W^p 



where VF(q) (or in the degenerate case) is an arbitrary function. A particular solution 

to the field equations can be retained only if the corresponding values of a and j3 satisfy 
(3). For a generic choice of W ^ some of the asymptotic AdS symmetries are broken. When 
the gravitational theory admits a field theory dual, the choice of boundary conditions (3) is 
mapped into a deformation of the action of the boundary theory [10, 11]. 

Adopting a convenient gauge, the most general spherically symmetric solution of the 
field equations can be written as: 

ds^ = -fir,t)e-''(^''Ut' + -^^+r'dn,.2, (4) 

f{r,t) 

If = f{r,t). (5) 
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The field equations stemming from (1) read: 



d,(^e'r'^)+^dr{r''-'e-'f^')-e-'V^ = 0, 

5' + ^ (r vv + v^'^) = 0, 

{d-3)il-f)-rf' + rf6'-^r^V = 0. 
In the static case, the field equations take the form: 

(^-3)(1-/)-./'-^/.'^ - S' 

/2 



6' 



np ■ 

T 



A constant shift in 6 corresponds to a rescaling of the time coordinate, therefore 6 appears 
only through its derivative in the equations for static solutions. For this reason, we imply 
the gauge choice: 

S{r) — > , 

for which gtt drr —1. Asymptotically, the field equations give: 

^+^^-^ f = ^ +l+ (rf_2),2A^-2 -;:d33+--- ^ 

2 Aa^og^r a(2A/3-Q)logr Mp 

+ + ^ + (d - 2)r'^-3 ^ (d - 2)r'i-3 r'^'-S + ' ' ' ' 

where Mq is an integration constant and a, f3 are the constants appearing in (2). The most 
general solution of the static field equations with 93 = is of the Schwarzschild-AdS (SAdS) 
type: 

/ = ^ + r^-^s^ Mo = rt' + rt-\ (6) 
5 = 0. 

However, for a variety of potentials V{(p) and spacetime dimensionalities, both analytical 
and numerical black hole solutions have been found for which the scalar field has a non- 
vanishing profile, and — t- at r — )■ 00. For a given potential, these solution generally form 
a two-parameter family spanned by the values (ph of the scalar field at the horizon and of 
the horizon radius: 

f{rh) = 0, ip{rh) = ifh . 

Alternatively, each solution is uniquely characterized by the (a, (3) coefficients of the scalar 
field, and the field equations provide a mapping {rh,iph) -H- {a, 13). Varying (ph and keeping 
rh fixed, one obtains a curve /3r^(a) (or ar^(/3)) in the (a,/3) plane, that is generally not 
single-valued. Finally, hairy solitons complete the ensemble of static, spherically symmetric 
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solutions to (1) with 7^ 0. These are defined by the requirement of regularity at r = 0, 
and form a one-parameter family spanned e.g. by the value of the scalar field at the origin: 

= ^{r = 0) . 

Correspondingly, the integration of the field equations provides a single /3s (a) curve (or 
as{/3)) describing the asymptotic parameters of the hairy solitons. In all the examples we 
have considered, one recovers /3s (a) (resp. as(/3)) as the r/j — )• limit of f3r^{a) (resp. 

The choice of boundary conditions at infinity restricts the ensemble of admissible hairy 
solutions to a one parameter family. For example, the number of admissible hairy black 
holes of radius is given by the number of intersections of the boundary conditions curve 
W^a{ct) with f3r^{a) or, in the degenerate case, of —W^p{l3) with ar^{(3). 

Asymptotic time translational symmetry implies the conservation of the total mass [12], 
whose expression receives a finite contribution from the slow falloff of the scalar field: 

A+/A_: M = fci(^il_^Mo + A_a/3 + (A+-A_)ty(a)) , (7) 

A+ = A_: M = ^(^Mo + ^/3^ + T^(/3)) . (8) 

The boundary conditions are insensitive to a constant shift in the value of W: the com- 
mon choice for boundary conditions that allow non-hairy solutions {W'{0) = 0) is to set 
W{0) = 0, in such a way that the mass of pure AdS vanishes. Formulas (7), (8) can be 
straightforwardly generalized to the case in which the boundary conditions curve can be put 
into the form (3) only piece-wisely: a suitable integration constant must be chosen on each 
piece, in such a way that W is continuous along the curve itself. 

As suggested in [13], for any one can build a constant entropy effective potential Vr^- 

MsAdsjrh) _^ A+ - A- 

A. = A_: VrM = + + f (10) 

where MsAdS is the mass of the SAdS solution of the same radius: 

(d-2)(r^-3+r^-i) 

Ms AdS = ys<i- ^ 



A+/A_: = '^^'^^+^^^(w{a)- ^Ja']^ (9) 



levrGrf 

Each extremum a* (resp. /3*) of Vr^ corresponds to a hairy black hole of radius satisfying 
the boundary conditions set by W . Moreover, the total mass of such black hole is simply 
given by Vgd-2Vry^{a*) (resp. Vgd-2Vrf^{fi*)). An effective potential Vs can be built out of 
the soliton curve in the same way: 

A+ / A_ : Vs(a) = {w{a) - £ ^ da'^ , (11) 

A+ = A-: Vs(/3) = ^^[w{p) + j\sdp'^ . (12) 

Vs plays a key role in giving the energy of the ground states: in [14], the authors have proven 
that the total energy in the theory (1) with boundary conditions (3) is bounded from below 
if and only if Vs has a global minimum. If this is the case, the minimum energy solution is 
precisely the hairy soliton that corresponds to the minimum of Vs • 
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Scalar, spherically symmetric perturbations of the "hair" 

We now turn to the possibility to extract additional information from the shape of the effec- 
tive potentials Vs and Vr^ • We start from a very simple idea, namely that some perturbative 
properties of the hairy solutions may be encoded in the form of the effective potential near 
its extrema. For this reason, we focus on the scalar perturbations of the hairy black holes 
and solitons. In [8], the authors have computed the frequency of growing modes of four 
and five dimensional unstable hairy black holes. In the remaining sections, we generalize 
and complete their results to models and boundary conditions that also admit perturba- 
tively stable black hole solutions, focusing on the dependence of such frequencies on the 
parameters that specify the boundary conditions. 

Due to the non-zero vev of the scalar field its perturbations give a backreaction on 
the metric already at the linear order. We limit our analysis to the case of spherically 
symmetric / = perturbations, that fit in the non-linear ansatz (4), (5): in this case, the 
linearized field equations can be reduced to a single master equation for the scalar field tp. 
Therefore, we consider the linearized ansatz: 

/ = /o(r)+eF(r)e-*-*, (13) 

S = 5o{r) + e A{r) e-'^^ . (15) 

Except for static perturbations, that describe infinitesimal displacements along the one- 
parameter family of hairy black hole solutions, the linearized Einstein equations let one 
express F and A' in terms of <I>: 

_2rfo^ 
d-2 ' 

A' = 

d-2 

The first of these equation can be used to extract the time-dependent variation of Mq as a 
function of the variations of a and f3, and check that the total mass (7), (8) is stationary 
around static solutions. Introducing the tortoise coordinate dr* = fQ^e^°dr, the master 
equation takes the familiar Schrodinger form: 

"^'^ foe-''Wo^ = -w^O) , (16) 



dr 



*2 



4r f{d-2){d-A) 3d - W ,2 ,4 



-/o [j^fo 



2r ,2 d-2 
2^ 



Near the horizon, located at r* = —00, /o vanishes exponentially, and the field variable 
satisfies a free equation: 



= e-*'^^ (1 + 0(6^'^'^ )) + <^out e*"^ (1 + 0{e'^' )) , (17) 

where k = fQ{rh)e~^'^ /2 is the surface gravity of the hairy black hole. Quasinormal frequen- 
cies are defined as the values of uj such that the master equation (16) admits a solution that 
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is purely ingoing at the horizon {^out = 0), and satisfies the boundary conditions at infinity, 
where the field perturbation decays as: 



A_ = A 



^ = Ar 



■+ ■ 



^ = Ar 



At the linearized level, the boundary conditions at the AdS boundary (3) read: 



A_/A+: B = W,aa{ao)A, 
A-=A+: A = -WMMB. 



(18) 
(19) 



Normal frequencies of hairy solitons are studied in the same way: ingoing boundary con- 
ditions at the horizon are replaced by the requirement of regularity at the origin. The 
corresponding Schrodinger problem is hermitian, and the eigenvalues need to be real. 
Therefore, in the absence of numerical issues, one finds only pure real or imaginary fre- 
quencies. In all the examples we have studied, as already observed for SAdS black holes 
[1] and many others, these normal frequencies coincide with the r/^ — )• limit of the quasi- 
normal frequencies characterizing the hairy black holes which correspond to the solitons 
(and compatible with the same boundary conditions). 

In the models we have considered, the background hairy black holes are known only 
numerically. Therefore, the generalization of the recurrence relation technique [1] inspired 
by the polynomial form of the SAdS metric (6) involves a recurrence containing infinite terms 
at each step. One needs to truncate not only the number of iterations, but the number of 
terms involved at each iteration [15], and the very interest of the original procedure is 
essentially lost. For this reason, our numerical technique consisted in a direct numerical 
resolution of the Schrodinger equation by shooting both from the horizon and the AdS 
boundary, imposing the relevant boundary conditions: ingoing propagation and (18), (19), 
respectively. Via a local optimization algorithm, one explores the complex plane of uj until 
the corresponding solution satisfies the correct boundary condition at the opposite end of the 
resolution interval ((18), (19) and ingoing propagation, respectively). Solving for uj in these 
two independent ways provides a solid consistency check to the numerical computations. 

As clearly explained in [16], the shooting method is severely limited in determining 
quasi-normal frequencies with imaginary part much larger than n in absolute value. This 
is typically the case for high overtones, and is due to the fact that the corresponding modes 
vanish rapidly near the horizon. Going back to the usual radial variable, near the horizon 
one has: 



The larger negative Im(a;), the more difficult is to resolve the presence of the outgoing mode 
near the horizon. To tackle this problem, the authors of [16] propose a resolution algorithm 
based on the relaxation method. Our work being oriented to the study of the parameter 
dependence, we have decided to implement the shooting algorithm with a few improvements: 

• When shooting from the horizon, the outgoing mode must be eliminated by imposing 
boundary conditions on the field variable ^ and its derivative This results in two 
distinct problems. First, depending on the value of Im(a;)/K, the outgoing mode can be 
subdominant with respect to the subdominant terms of the ingoing mode. Therefore, 
we compute analytically, by deriving the master equation, a sufficiently large number 



<^ = ^in{r - r;,)— (1 + 0{r - m)) + ^out{r - r^r'^'' (1 + 0{t - r,,)) . 
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of coefficients of such subdominant terms. To be precise, one must know the value of 
n such coefficients, where: 

Still, the numerical precision must be sufficient to resolve all such subdominant terms 
in the vicinity of the horizon. This turns to be the decisive constraint in the compu- 
tation of high overtones. 

• When shooting from the boundary of AdS, one can fit the solution near the horizon 
in order to estimate ^out and <l>j„. This method has the advantage not to require 
the evaluation of the coefficients of the subdominant terms of the ingoing mode: one 
can fit, besides than ^out and <^>j„j the coefficients of a the power expansion near the 
horizon: 

$ = ^^n{r-rh)-^'^l^''{l + Cl{r-rh)+C2{r-rhf + ... + Cn{r-rhT) 

However, the limitation imposed by the numerical precision of the procedure is always 
present. 

Depending on the model and the boundary conditions, it has been possible to find the first 
two or three overtones with n = 1 or n = 2. 

The repeated procedure consisting of solving the field equations, find the hairy solu- 
tions compatible with certain boundary conditions and then computing the quasi-normal 
frequencies that characterize them, has been implemented as a Mathematica library con- 
taining simple optimization algorithms and methods to manipulate the input/output. The 
library is adapted for general spacetime dimensionalities and models of the type (f), and is 
available from the author upon request. 

3 Hairy solutions in D = 5, A/^ = 8 gauged supergravity 

First, let us consider the following model: 

where 7 = y^2/15, that emerges as consistent truncation of five-dimensional J\f = 8 gauged 
supergravity [5, 17]. In this case, the Breitnerlohner-Freedman bound = —4 is satu- 
rated: 

a logr /3 

Gauged supegravity is thought to be a consistent truncation of ten-dimensional type IIB 
supergravity so, with usual boundary conditions a = 0, we know the dual theory to be 
M = 4 Super- Yang-Mills on M x 5^. The bulk scalar ip couples to a trace operator of the 
scalars in the adjoint: 

I3 = {0), = cTr U?-^ J^c&n . (21) 
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Fi gure 1: a,-^(/3) curves describing hairy black holes of different radii. In the right panel, log-log 
plot of the /3 > 0, a > quadrant: all the a,.^ converge to the asymptotic Uc value with a seemingly 
r/j -independent power law for large r^. The gray zone contains no black hole solutions. 



Imposing generalized a = —W^p boundary conditions corresponds to deforming the dual 
SYM action [10]: 

SsYM ^ SsYM - ^ j dt d^9 Vf W{0) , (22) 

where F is the metric on unit S^. 

For positive ip, the bulk potential V is concave and one can find numerically a hairy black 
hole solution for any value of and (ph > 0, with the scalar field vanishing at infinity. At 
if = if^ < 0, the potential has a local minimum that corresponds to a smaller radius AdS, 
and hairy black holes can be found only for ip^ > ip^ : below this value, the potential starts 
increasing exponentially, and the the resolution of the field equations leads to (/J — )• — oo at 
a finite value of the radial coordinate. Some curves are plotted in Figure 1: for positive 
large iph, ctr^ converges to an asymptotic value ac — 0.54772, via a power law behavior 
that seems to be radius-independent at large rh (right panel). On the contrary, (3 becomes 
arbitrarily large. In this region, no black hole is found with a < ac- When <C 1, as we 
have pointed out in the previous section, a^^ coincides with the as curve described by hairy 
soliton solutions. At large (3, as converges to ac with a logarithmic error that is small for 
/3>10. 

For negative values of (ph, one finds /3 > and a < 0. When approaches (p- from 
above, both /3 and a diverge with a linear behavior corrected by logarithmic terms in /3: it 
takes more and more "time" (radial coordinate) for the scalar field to escape the fixed point 
represented by the smaller radius AdS solution, that corresponds to constant negative ip, 
thus to infinite negative a. However, following previous studies, our analysis has focused on 
the hairy solutions with ip > 0. 

We shall consider the following boundary conditions: 

W = -{/32 + ^/3^ 
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Figure 2: In the left panel, intersection of the boundary conditions curve a/^^(/3) (/ = 10~^, 
77 = 0.2, red daslied curve) with some ar^{P) curves. Each intersection point corresponds to a single 
hairy solution compatible with the boundary conditions. For every < rf^, one finds two distinct 
hairy black holes: the two resulting branches are plotted in thick, solid lines. In the right panel, plot 
of the corresponding constant entropy effective potentials (SttGs = 1), whose extrema correspond to 
the same hairy black hole solutions: at r/j = r'f^, the two non-trivial extrema disappear. 



where / and rj are positive parameters. Similar boundary conditions were considered, in a 
very different perspective, in [9]: the dynamics at short times of the dual field theory was 
explored in the limit 77 —t- 0, in which suitable initial conditions evolve into a cosmological 
singularity. When rj is finite, it was argued that the same evolution ends with the formation 
of a large hairy black hole [5] , of the kind we are going to study. 

When / is small, say / < 0.1, a/^^ crosses the soliton curve at /3 > 10 where Os — cuc- 
Hence, the effective potential for solitons (12) takes the approximate form: 



When T] is small enough: 



^^^e=^-0.495 



Vs has two non trivial extrema fi± that correspond to hairy soliton solutions compatible 
with the required boundary conditions (see Figure 2). The less hairy soliton, located at local 
maximum /3_ of V^, has positive mass. The more hairy one, located at the bottom of \>s 
has negative mass provided 77 < f?c/2- Two distinct branches of hairy black holes start from 
these solitons, and end by joining to each other at some finite radius = r^- Accordingly, 
while r/i grows, the shape constant entropy effective potential changes continuously and, 
for r/i > no more extrema are present. In between this critical point and the solitons, 
above some intermediate value r™'^, the more hairy black hole becomes more massive than 
the Schwarzschild-AdS (SAdS) solution of the same radius. 
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Quasi— normal frequencies 

The shape of the soHton effective potential suggests that the less hairy soliton may be 
unstable under perturbations of the scalar field. Such instability has indeed been found in 
[7] with simpler boundary conditions: 

a = a/,0 = fP , 

that give rise to a soliton potential that is unbounded from below. The authors have found 
that the corresponding black holes possess a growing mode with pure imaginary frequency 
u = iuQ, I'o > (see (13)). 

One can actually prove that, if such instability exists, the associated frequency ought to 
be pure imaginary. To simplify the notation, let us choose r*(r = oo) = 0. Multiplying (16) 
by and integrating between r* = — oo, corresponding to the horizon, and r* = — e, one 
finds: 

' dr* ('^*f^ - foe~''"Vom' + oo^mA = . 

For a growing mode of frequency uj = ujr + iuji, ujj > 0, ingoing boundary conditions at the 
horizon: 



imply that <5 decreases exponentially when r* — t- — oo, so the integral above is perfectly 
finite. Taking the imaginary part and integrating by parts, we obtain: 

Im ( ^*^J (r* = -e) + 2ujruji J dr* = . 
The coefficient function foe~'^^°Vo has, near the AdS boundary, the following behavior: 

foe-^''Vo = -^\r*\-^ + Oil) . 

Hence, asymptotically: 

$ = -^|r*|^/2 log |r*| +S|r*|^/2 flog |r*| \r*f/'^ 



Imposing our boundary conditions (19) requires A and B to have the same phase, so near 
the boundary: 

Im(<&*^)=0(log|r*||r*p) . (23) 

So, moving to the limit e — )• 0^ we deduce ujr = 0. 

Going back to our af^r] boundary conditions, we observe that when 77 is small, the shape 
of the effective potential around its unstable maximum is essentially r/-independent, so 
we expect to find the same instability: our numerical data confirms that a growing mode 
characterizes the less hairy branch of black holes with aj^^ boundary conditions as well. 
This frequency is accompanied by another pure imaginary but damped frequency uj = wq, 
uq < 0. In particular, we find ^0(^^/1 = 0) = —1^0(^/1 = 0), and this comes as no surprise: 
a given regular solution of (16) solves both for uj and —uj, translating the symmetry of 
both the wave equation and the boundary conditions under time reversal. At finite radius 
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0.0 0.5 1.0 1.5 2.0 0.0 0.5 1.0 1.5 2.0 



n n 

Figure 3: f — 0.1, 77 = 0.2. Real and imaginary part of wq, wq, ivo (solid lines, same colors as in 
Figure 2) and further quasi-normal frequencies uji and UJ2 (dotted lines). Positive imaginary part is 
associated with growing modes. At the critical point rh — f^, and ivQ vanish. For a more detailed 
plot of the critical point r^, = r^, see Figure 4. 

r/i > 0, an ingoing solution still solves for —u too, but it is purely outgoing with respect 
to this frequency^. Accordingly, we find vq —I'o at finite radius. 

When the radius is increased and approaches the critical value r^, fo('^/i) approaches zero 
(see Figure 3). Indeed, at r = r^, the boundary condition curve aj^ri and 0,.= are tangent to 
each other, and the static perturbation giving rise to the infinitesimal static displacement 
along this tangent is compatible with the boundary conditions. On the contrary, the spurious 
frequency ii'oirh) stays away from zero at r = rf^. These frequencies can be followed back 
along the branch of more hairy black holes, that starts at r/j = from the soliton that 
corresponds to the global minimum of V^. We will denote ujQ{rh) and ujQ(rh) the related 
frequencies (see Figure 4): 

Mrl) = iMrh) = , 

For values of the radius extremely close but smaller than r^, loq and are still pure 
imaginary, but uq has now a negative imaginary part corresponding to damped oscillations. 
Within a typically tiny distance from r^, ujq joins with a)o (see Figure 3) at some point of 
the imaginary axis. For smaller values of rh, the remaining frequency acquires an oscillating 
part Re(6t;o) ^ and is now distinct from his partner at uj = —ujq- Finally it ends, at = 0, 
in a pure real normal frequency of the more hairy soliton. A very similar bifurcation has 
been observed in the analysis of vector quasi-normal frequencies of D7 branes [16], with 
momentum replacing as driving parameter. In that context, the presence of a = 
mode at A: = translates the homogenous diffusion of a conserved charge in the presence of 
translational symmetry. 

^On the contrary, is an ingoing solution for —uj*, and the quasi-normal frequencies come as usual in 
couples (tj, —LJ*). 
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Figure 4: In the left panel, trajectory of ujq, ujq, ivq and H^q in the complex plane; arrows correspond 
to the direction of increasing radius, and the blue points mark the critical radius where wq = 
Uq = 0. In the right panel, a detailed plot of Im(wo) and Im(wo) near r^, = rj^; Re(cjo) is plotted in 
dotted lines. At the bifurcation point, the ujq frequency characterizing the stable branch becomes 
pure imaginary and gives rise to two distinct pure imaginary frequencies loq and cjo, one of which 
turns to zero frequency when — >■ rf . 



In summary, with Oif^jj boundary conditions, the qualitative picture given by the evo- 
lution of the constant radius effective potential (see Figure 2) turns out to be accurate in 
predicting the perturbative stability/instability of the hairy black holes under perturbations 
of the hair itself. 

Apart from the one we have just described, both branches of black holes possess a 
presumably infinite spectrum of other quasi-normal frequencies. They all correspond to 
damped perturbations and, as we have already pointed out, they are characterized by pure 
real frequencies in the limit — t- 0. By continuity between the stable and unstable branches, 
their trajectory as a function of the radius forms closed curves in the complex plane (see 
Figure 3). In the range of parameters we have explored, the sequence of these "standard" 
quasi-normal frequencies respects the order one finds, for example, for standard perturba- 
tions of SAdS black holes: they can be labeled in such a way that both Re(a;„) and |Im(a;„)| 
increase with n and this, for any r^. On the stable branch, ujq does not always respect such 
order. In fact turns out to be the longest lived frequency but, for small values of r], its 
real part can be greater than some of the Re(a;n). 



Parameter dependence and effective potentials 

Let's consider again the approximate expression of the soliton effective potential in the 
/3 > 10 region: 



where we have made use of the standard AdS/CFT relation: 
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Recalling the form of the dual operator (21), we observe that the linear term in Vs(/?) 
takes the form of a mass term for the scalars. Such kind of term is already present in the 
undeformed SYM action as conformal coupling of the scalars to the curvature of : 



SsYM= / dtd-'eVT 



2 2 V SvrGi 



+ .... (25) 



Moreover, a classical analysis shows that the extrema of the SYM deformed scalar potential, 
that may be identified with our hairy solitons, are characterized by <I'i = for i = 2, . . . , 6. 
Therefore, comparing (24), (25) and (21), we tentatively identify: 



27r2 



a, 



.iV2 



(26) 



At this stage, this is a simple prescription for the soliton effective potential to coincide with 
the deformation plus a free term in the original SYM action. The masses of the hairy solitons 
can now be obtained as vev of the SYM scalar potential, integrated over the volume of S^. 
However, such prescription can be checked by an independent calculation, e.g. matching 
the normalization of the correlation function of O in the field theory. As we shall show in 
the next section, the mere existence of such prescription is not guaranteed, but is rather a 
property of the particular consistent truncation we are considering: in fact, the key point 
is that as — >• Oc when /? — t- oo. It should also be stressed that, as the authors of [7] have 
pointed out in a four dimensional case, the mass term does not appear at small /3, where Og 
grows as 

We can cross-check our prescription via the quasi-normal frequencies data. We introduce 
a further simplification by tentatively neglecting the non-abelian structure of the Higgs 
scalar <I>i: 



$1 = VNcI)In , (27) 

P = —cP'. (28) 
ac 

Neglecting all the other fields in (25), we consider the following toy model free action for 
the dual theory: 

S = N'Jdtd'9Vf(^-l{d,ct>r-l^'-^^W,,,(^^4>')) (29) 

= N'Jdtd'eVf ^^-lidc^f-^^ . 



At the extrema of V^, the mass of the scalar field (f> is given by: 

.2^,„ad^^s _ 2/3 fda 

In particular: 



"^'^4c/3^ = ^(^ + T^,^;3) . (30) 



m2(/3+) = t^2>o. 

We propose to compare these masses to the values of —Vq and loq that characterize the 
unstable/stable solitons respectively. We generalize this proposal to match the values of 
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Figure 5: In the left panel, imaginary part of uq as a function of the radius, for ?/ = 0.2 and 
/ = 10^^,10^^ and 10^'^. The cusps coincide with the bifurcation point (see Figure 4). When 
/ — 0, corresponding to the large field limit, the ujq becomes stable. In the right panel, the goal 
function whose roots coincide with the quasi-normal frequencies, for the hairy black hole of radius 
rji ~ 2.4 and compatible with boundary conditions r] — 0.2, / = 10"^: the ujq frequency has an 
imaginary part of order — 10~^, while the second overtone uji has Im(a;i) ~ —12. 



and lOq for > 0: we already know that, when r/j approaches r£, one finds 1^0,000 — )• 0. 
At the same time, the cmvatm'e of the effective potential Vr^ near its unstable maximum 
vanishes (see Figure 2). This suggests the minimal generalization of (30) that consists in 
replacing Vg with V, 



rh- 



2 , r.d'^Vr, 2/3 /dar, 



and accordingly: 



The derivative dar^/df3 can be easily evaluated numerically. Of course, we cannot hope to 
describe in this simple way the dissipation described by the imaginary part of ujq on the 
stable branch of hairy black holes: is real while ujq is not. Therefore, we limit ourselves 
to comparing u;^^ with the real part of Wq as a function of the horizon radius. 

The numerical agreement depends on the values of / and r], and it can be better analyzed 
by taking a closer look at the parameter dependence of the toy model soliton potential for 
(p that appears in (29): 

where / and fj include the order one normalization factors in (29). Of course, is simply 
given by the second derivative of V^^.^ with respect to <j), evaluated at his extrema. The way 
we parametrized the boundary conditions lets one reabsorb / in a rescaling of and V/^^: 

VfJ(t>) = r'ViJy^<P). (31) 
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Figure 6: Comparison between the real part of Wq and the curvature of the constant entropy 
effective potential at its extrema, for boundary conditions with t] = 0.2 and / = 10~^ (top left 
panel), / = 10~^ (top right), / = 10^^ (bottom left), / = 10"" (bottom right). Scattered points 
are quasi-normal frequencies (stable branch) and Wq, while green dots represents the damped, 
pure imaginary frequencies zwq characterizing the unstable black holes. Dashed lines represent the 
second derivative of the radius-dependent constant entropy effective potential, evaluated at its two 
extrema. In the large field limit / — >■ 0, the agreement becomes better and better, and = i/q. 
Similar results have been obtained for smaller values of rj. 

Thus, when / decreases, the masses of the stable/unstable solitons and the corresponding 
values of cf^ oc (3± scale as f~^. At the same time, as one can deduce from Figure 1, the 
critical radius increases; numerical data shows that: 

rl oc . (32) 
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From (31), one sees that uj and do not scale with /: 



Therefore, as long as this correspondence holds, a;g(r/i = 0) and i^o(''^ ~ 0) should reach 
a finite limit when / — )• 0. This is confirmed by numerical results (see Figure 6), and also 
holds for hairy black holes, provided one rescales according to (32). Moreover, in the 
same limit / — t- 0, we observe that that the imaginary part of ujq rapidly goes to zero along 
the whole branch of stable hairy black holes (see figure 5). One finds |Im(a;o)| oc f^'^ at 
constant at constant radius, |Im(ti;o)| oc at constant temperature. This justifies the free 
character of our toy model action, so as to study the cuq frequency. In a related way, the 
spurious frequency ii)o characterizing the unstable black holes appears to coincide with the 
opposite of its companion ivo, signaling again the approach to a non dissipative regime: 

Finally, the values of the curvature of the constant effective potential near its extrema, 
and —i^r^, give in this limit a very good estimation of Re(a;Q) and — fg (see Figure 6). 

In summary, / — t- corresponds to increasingly large values of the condensate f3. This 
limit preserves the curvature of the deformation in terms of boundary scalar (p: the masses 
of the solitons scale as f~^, but so does (/>^. We have reviewed the presence of a particular 
frequency ii'q (resp. uq), in the spectrum of quasi-normal frequencies of the corresponding 
unstable (resp. stable) hairy black holes, and showed that its value is very accurately 
predicted, in the limit / — t- 0, by the value of the curvature of the constant entropy effective 
potential at its maximum (resp. minimum). More specifically, ii^o is the frequency of the 
growing mode that is responsible for an instability of the related hairy black holes. This 
agreement is realized by assuming that the mass term for /3 appearing in is the conformal 
mass term coupling <I>i to the curvature of S^: under this hypothesis, the normalization of 
C, that enters in the expression of rn^^, is fixed by the numerical value of Oc- In the same 
/ — )• limit, the spurious frequency iDq approaches the value —ivq: therefore, these two 
frequencies appear as solutions to the toy model equation for homogeneous perturbations: 



iV2 dd)-^ 



0, 



that generalizes (29) to finite radius. At the same time, the imaginary part of ojo goes to 
zero, and ujq solves the analogous equation: 



I 1 d^Vr, 



</>+ 



While Wo is longer and longer lived when / decreases, the remaining "standard" quasi- 
normal frequencies do not share this feature: for example, as we show in figure 5, the 
imaginary part of toi reaches larger and larger negative values, as one usually expects for 
increasing values of r^. 

In conclusion, the normalization (26) is quantitatively confirmed by the agreement of 
the quasi-normal frequencies uoq^ivq with the curvature of the constant radius effective 
potentials, in the limit of large condensates / — )• 0. Therefore, we expect — )• Oc 7^ 
to be a consequence of the curvature of the spatial sections of the boundary. In the dual 
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field theory, homogenous excitations of the condensate operator around the minimum of the 
effective potential are longer and longer lived, due to the presence of a single quasi-stable 
resonance. Remarkably, higher overtones do not share this property: their imaginary part 
remains finite when / — )• 0. 

4 Hairy solutions in D = 4, J\f = 8 gauged supergravity 

We now turn to analyze four dimensional hairy black holes, starting from a generalized 
consistent truncation of four-dimensional of = 8 four dimensional supergravity [14, 18]: 

- /^'^^ (f 4<^^>= (^) (^) (^)) . 

where b and c are free parameters and: 

n = \/l + 62 + c2 . 

The terms in the potential are combined in such a way that m? = —2. Accordingly, the 
scalar field decays at the AdS boundary as: 

a p 

and we consider boundary conditions of the form /3 — ^^^a • 

With usual boundary conditions a = 0, the dual field theory is defined on M x S^, and 
it has been conjectured to be the so-called ABJM theory [19]. a now represents the vev of 
a dimension one operator, and boundary conditions /3 = W^a induce a deformation of the 
boundary action: 

a = {0), s^s--^ dtd^eVfw{o). 

87rG4 J 

Like in five-dimensions, one can obtain hairy black hole and hairy soliton solutions by 
numerically solving the static, spherically symmetric field equations. The ensemble of solu- 
tions of a given radius r/j determines a curve in the (a, (3) plane, that we now parametrize 
as /3j.^(a). The scalar potential being even, we can restrict our analysis to solutions with 
Q > 0. As usual, when r/^ — t- one recovers the curve characterizing hairy solitons solutions, 
/3s (a). Solving the scalar field equation on AdS one can show that: 

/3^(a = 0) = -- . 

TT 

At the non-linear level, the precise form of /3s clearly depends on the particular truncation. 
For generic values of b and c the soliton curve scales, at large a, as one would expect from 
dimensional analysis [14]: 

Ps^-ScO^, Sc>0. 

According to the definition of Vs (11), this implies that the soliton effective potential contains 
a positive definite term proportional to a^. This term insures a lower bound to the energy 
even in presence of boundary conditions that imply a perturbative instability of AdS [14, 20]: 

W{a) = -{a\ />- 

Z TT 
/?/ = -/«• 
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Figure 7: In the top panels, some (o^) curves for 6 = c = and 6 = c = 1 respectively; the 
black, dotted line represents the slope of /3s at the origin, independent from b and c; the red, dashed 
lines represent /?/ boundary conditions for / = 0.25 and / = 1 respectively. In the bottom panels, 
we show the corresponding form of the constant radius effective potentials. For b — c = 0, the hairy 
black holes corresponding to the maximum of Vr^ become more and more hairy and massive, with 
respect to SAdS, as grows. For b — c — 1, the hairy black holes have lower mass than SAdS, until 
they disappear at — rf^- 



When / is large, the soUton effective potential takes the mexican hat form: 

The hairy soliton sitting at the minimum of Vg has negative mass. When the radius is 
increased, the slope of the origin becomes more and more negative and eventually 

P'r^ioi = 0) < — / (see Figure 7, right column). Hence, one finds hairy black hole solutions 
compatible with f3f only for rh smaller than some critical value r£ or, equivalently, for 
temperatures smaller than some T^. Both and are increasing functions of /, that can 



19 



be evaluated analytically for large /: 



(33) 



In the picture provided by constant radius effective potentials, the non trivial minimum 
of Vr^ disappears when approaches rj^ (see Figure 7). In the minimal generalization of 
the bulk action containing a f/(l) gauge field, this is interpreted in the dual field theory as 
the continuous restoring, at high temperature, of a spontaneously broken symmetry. The 
breaking is induced by the deformation corresponding to the boundary conditions, without 
the need to add a finite chemical potential [6]. 

The truncation obtained with 5 = c = provides a completely different picture. As 
already noticed in [7], here the soliton curve approaches a constant value when a goes to 
infinity, /3 ^ f3c — —0.7071. This asymptote corresponds to a mass term for O in the soliton 
effective potential. Adopting Pf boundary conditions, the latter takes the asymptotic form: 



Therefore, one finds a hairy soliton compatible with /3j boundary conditions represented 
by the maximum of Vs (see Figure 7, left column), that is now unbounded from below. At 
finite radius, the /S^^ curve converges to a constant value that, unlike in five dimensions, 
increases with r/j. Therefore, f3f boundary conditions admit hairy black holes of arbitrarily 
large radius. 

For definiteness, we shall concentrate on the models defined by 6 = c = 1 and 6 = c = 0, 
adopting 13 j boundary conditions. 

Truncation b = c = 

For 6 = c = 0, the shape of the effective potential suggests again the presence of an 
instability of the hairy solutions under perturbations of the hair. Such instability has indeed 
been found [8] with different but qualitatively similar boundary conditions, and appears a 
single pure imaginary frequency with positive imaginary part in the spectrum of normal 
perturbations. Like in five dimensions, one can again show that the frequency of such a 
mode must necessarily be pure imaginary. The proof follows the same line as before: near 
the boundary, the reduced perturbation variable behaves as follows: 



where now the boundary conditions impose B = W^aaict)A. The last step of the proof (23) 
is replaced by: 



The presence of this perturbative instability is confirmed by our numerical results: an pure 
imaginary frequency uj = ivq, > Q characterizes all the hairy solitons and black holes 
we have considered. Again, a spurious frequency Cj = ii>o, i^o < is also present, with 
t^oifh = 0) = i'Q{rh = 0). When / goes to zero, both fo(r/i) and i^oirh) converge to a single 
universal function of the radius (see Figure 8)'^. This suggests that, like in five dimensions, 

^This feature also appears with the j4d5-invariant boundary conditions /? = — jo? considered in [8]. This 
has probably not been recognized by the authors due to a background dependent gauge choice, &{rh) = 0: 
this implies that 5 has different values at the AdS boundary for different backgrounds, where times are 
measured in different units. 




(34) 



^> = yl-Sr*+0(r*2) , 




20 







2 



4 6 



8 



Figure 8: b=c=0: Comparison between the squared frequency —Vq of the growing mode and 
the curvature of the constant entropy effective potential, Vj-^^, for various values of /: from the 
outermost curves, / = 0.5, 0.25, 0.125, 0.0625, 0.03125. The red dashed line at the center represents 
the curvature of the effective potential Vr^ . In the large field limit / — 0, one finds Vq=Vq= ~^rh ■ 

in the limit of large condensates the frequency of these perturbations may be described by 
the curvature of the constant radius effective potential. Using the approximate expression 
(34), we see that the soliton solution is located at a_ oc /~^. Therefore, considering a 
boundary scalar (p such that a oc 0^, we find again: 



However, to make a quantitative comparison, we need to fix the normalization of a in 
terms of (p. Like in five dimension, we tentatively focus on the mass term appearing in the 
expression of the soliton potential at large a: 



As a represents the vev of a dimension one operator O, and the dual field theory is defined 
on M X S"^, we try to identify such term with the conformal coupling of O to the curvature 
of S'^ [7] . Neglecting the non-abelian structure of the field theory, we again consider a toy 
model free theory for a single scalar field: 



87rGiVs ^ \(3c\a + W{a) . 



O 



\Pc\a + W{a) 



1 




W{a) 



s 



87rG4 

7rG4 
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c 



|/3c 
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where 4> is supposed to be canonicahy normaUzed in the dual field theory. With this pre- 
scription, we can compare the quasi-normal frequency l^Q , clS 3< function of rh, with the 
curvature of the constant radius effective potential Vr^ at its extrema, computed in (f) vari- 
ables: 

,2 _ f^^K,, _ . _ d^Vr^ _ ( dPrn 



Once again, can be easily evaluated numerically, once the hairy black holes are found. 



In Figure 8, we show the numerical results for the agreement between and —Vq, —i^o- 



„U™ ^1, ;„„i „„„„U„ „„ ^ U„^,„ .^^^ ^^^^ 

While / decreases, the agreement gets better and better, and i^o ~ z>o as one would expect 
in the absence of dissipation. Therefore, like in five dimensions, this truncation exhibits 
a large field regime in which the effective soliton potential can be consistently identified 
with the multi trace deformation plus a mass term that couples the dual operator with the 
curvature of S'^. We have also considered "regularized" boundary conditions that recall the 
one adopted in five dimensions [13] : 

= -/a + /\a^ r/>0. 

The picture we have obtained is essentially the same: the regularizing term produces a new 
branch of perturbatively stable black holes that have negative mass at small radius. They 
are characterized by a quasi-normal frequency uq that, in the small / limit, becomes stable 
(Im(a;o) — ^ 0) and whose real part is well predicted by the curvature of Vr^- 



Truncation b = c = 1 

For 6 = c = 1, the shape of the effective potential with /3j boundary conditions strongly 
recalls the one we encountered in the study five dimensional hairy black holes. There, the 
critical radius corresponded to the joining point of the two extrema of Vr^ ) both located at 
/3 > and corresponding to unstable/stable hairy black holes. Here, SAdS black holes play 
the role of the unstable hairy black holes. Therefore, we expect to find a special frequency 
Wo, in the spectrum of quasi-normal perturbations, that vanishes for — t- r^. However, the 
situation is rather different from the previous cases. The large / form of V^: 

87rG4V, --^a^^ya^ , 

is dominated, when a is large, by the term ScO^, that is not present in the deformation 
VF(a) and cannot be identified as a free term of the undeformed field theory action. In fact, 
the behavior of the quasi-normal frequency shows a qualitative resemblance with the five 
dimensional case (for comparison. Figures 6 and 9). However, the limit of large condensates 
is very different: when / grows, one recovers the scaling suggested by dimensional analysis: 

wo oc /"^ , (35) 

both for the real and imaginary part of uq. In particular, cjq reaches the imaginary axis at 
a value of that scales as and is significantly smaller than r^. Therefore, even if the 
scaling (35) naturally characterizes the curvature of Vs for the same dimensional reason, we 
cannot expect to find a quantitative agreement between Re(a;Q) and . 
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Figure 9: b = c = 1: real and imaginary part of the least damped quasi-normal frequency of 
hairy black holes, for boundary conditions with / = 1 (top left panel), f — A (top right), / = 16 
(bottom left), / = 64 (bottom right). The imaginary part is also shown. Solid lines represent data 
from quasi-normal frequencies, while dashed lines represent the square root of the second derivative 
of the radius-dependent constant entropy effective potential, evaluated at its two extrema. The 
frequencies at small radius, relatively to /, turned out to be above the confidence limit of our 
numerical computation (see discussion in Section 2), and are not show. In the large field limit 
/ — > oo the real and imaginary part of cjq, as functions of f^^rji, globally scale as /. Therefore, ujq 
is still associated with a large dissipation. 

Conclusions 

We have given detailed numerical results for the I = scalar quasi-normal frequencies of four 
and five dimensional AdS black holes with scalar hair and generalized boundary conditions. 
The computation has been performed by taking in account the full backreaction of the scalar 
field on the metric. 

We have identified models in which, at large values of the "hair" , the difference between 
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the effective potential Vg and the deformation of the dual theory can be naturally identified 
as a conformal mass term in the field theory action. As a consequence, the normalization 
of the operator with respect to a canonically normalized scalar can be numerically obtained 
from the soliton data. The corresponding hairy black holes possess a quasi-normal frequency 
ojQ which is quasi stable (|Im(a;o)/'^o| 1) and coincides with the curvature of the effective 
potential at fixed entropy Vr^, provided the condensate is large enough. In the dual field 
theory, the corresponding mode is associated to homogenous perturbations of the condensate 
that undergo very weak dissipation. A heuristic interpretation for the presence such mode is 
that higher derivatives of Vr^, corresponding to self-couplings of the perturbations, vanish 
in the limit of large condensates. Oppositely, the frequency of higher overtones mantain a 
finite imaginary part. 

Boundary conditions have been considered for which two branches of stable and unstable 
hairy black holes join at a critical point, where the curvature of Vr^ vanishes. At the critical 
point, the mode mentioned above has zero frequency and describes the transition between 
the two branches. For smaller values of the radius, its frequency turns to pure imaginary 
values before bifurcating into two frequencies with non zero real part. 

Finally, we have discussed a four dimensional model in which Vg contains an interaction 
term deduced from the supergravity data, which is not present in the deformation. Here, 
the imaginary part of the least damped frequency does not vanish in the limit of large 
condensates, and the constant entropy effective potential describes its real part only at the 
level of the scaling with the parameters of the deformation. This quantitatively different 
behavior may be interpreted by the fact that, unlike in the first class of models, some higher 
derivatives of do not vanish in the limit of large condensates. 
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